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ABSTRACT :  The  main  purpose  of  this  work  is  to  perform  a  new  denoising  method  based  on 
a  nonlinear  anisotropic  diffusion  for  the  reducing  of  the  multiplicative  speckle  in  high 
resolution  Synthetic  Aperture  Radar  (SAR)  images.  In  order  to  be  applicable  to  the  sampled 
data  and  reduce  computing  complexity,  an  efficient  discretization  scheme,  e.g.  additive 
operator  splitting(AOS)  scheme  is  chosen  here  for  the  proposed  approach.  Examples  are  given 
comparing  the  new  filter  to  established  speckle  noise  reduction  methods  on  an  ERS-2  SAR 
image.  Experimental  results  show  that  this  method  not  only  can  effectively  reduce  speckle, 
while  preserving  important  discontinuities  present  in  the  input  data,  but  also  obtain  a  better 
performance  than  the  traditional  local  statistical  filter. 

1  INTRODUCTION 

The  basic  goal  of  SAR  images  denoising  is  that  preserving  the  image  detail  characteristics 
like  edges  and  textures,  while  suppressing  the  speckle  of  uniform  area  in  the  image.  The  most 
classical  approach  can  be  implemented  by  multi-look  despeckling  in  the  imaging  processing. 
But  this  increase  in  radiometric  resolution  is  gained  at  the  expense  of  spatial  resolution.  As  for 
the  multiplicative  characteristic  of  the  speckle,  many  spatial  domain  filtering  algorithm[l]  have 
been  developed,  and  gradually  become  the  mainstream  technique  of  the  SAR  image 
despeckling,  the  typical  approaches  include  the  adaptive  local  statistical  filter  (like  lee  filter, 
kuan  filter,  etc)  and  geometric  filter  etc,  the  key  of  this  kind  of  approach  consist  in  how  to 
achieve  effective  despeckling  according  to  the  statistical  property  of  the  sliding  estimated 
window  and  choose  the  appropriate  window  size  to  gain  the  biggest  homogeneous  region.  In 
addition,  filtering  technique  based  on  iteration  and  wavelet  approach  [8]  and  EMD  approach  [2] 
proposed  by  Han  et  al.  both  obtain  a  good  performance  in  reducing  speckle. 

In  recent  years,  the  image  processing  approach  based  on  nonlinear  diffusion  have  already 
been  applied  in  the  field  of  image  enhancement,  image  segmentation  and  scale  space  analysis 
etc  [3]  [4]  [5] [6],  its  trait  is  reducing  the  noise  and  sharpening  edges  at  the  same  time.  In  order 
that  the  continuous  anisotropic  diffusion  equation  is  applied  in  the  discrete  image  data,  the 
nonlinear  diffusion  filtering  algorithm  based  on  additive  operator  splitting(AOS)  scheme  is 
proposed,  and  then  compare  its  performance  with  the  traditional  SAR  speckle  filter,  the 
experiment  results  show  that,  based  on  several  evaluation  measure,  this  approach  not  only 
achieve  good  balance  between  reducing  noise  and  preserving  edges,  but  also  obtain  a  better 
performance  than  the  traditional  local  statistical  filter. 
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2  NONLINEAR  DIFFUSION  FILTERING 
2,1  The  Physical  Background  of  Diffusion 

Firstly,  giving  the  heat  diffusion  equation 

dtu  =  div(D  -Vu)  (1) 

Where,  u  is  concentration,  div  is  the  divergence  operator  ,V  is  the  gradient  operator  ,  D  is 
diffusion  tensor  and  t  is  time.  This  equation  appears  in  many  physical  transport  processes.  In 
image  processing,  the  concentration  can  be  considered  as  the  gray  value  at  some  position.  If  D 
itself  depends  on  the  evolving  image  u,  the  resulting  equation  describes  a  nonlinear  diffusion 
filter. 

Nonlinear  diffusion  was  first  proposed  by  Perona  and  Malik  in  1990[7],  and  built  the 
following  isotropic  diffusion  model: 

dtu  =  div(g(\  Vu  |)  •  Vu)  (2) 

But  Perona-Malik  model  still  have  some  fault.  Thus,  Catte  et  al.  proposed  that,  by 
substituting  the  gradient  Vu  with  its  estimate  \VGa*u\  where  Ga  can  be  any  smoothing 
kernel,  nonlinear  diffusion  filter  can  be  more  robust  for  noise.  While  as  for  the  problem  of  the 
unfiltered  edge  noise,  it  have  to  turn  to  the  anisotropic  filters  where  the  flux  j  =  -D  ■  Vu  is  not 
generally  parallel  to  the  image  gradient. 


2.2  Anisotropic  Diffusion  Model 


Generally,  the  presence  of  the  noise  near  the  edges  result  in  the  fluctuation  of  the  image 
gradient,  consequently  cause  the  fluctuation  of  the  gradient  intensity  and  direction.  Nonlinear 

isotropic  diffusion  filter  includes  a  scalar  diffusion  tensor  D  =  g(|  Vun  |) ,  so  the  flux  always 

parallel  to  the  image  gradient  Vu  ,  noise  near  the  edges  can  not  be  smoothed,  however  rotate  it, 
then  diffusion  filter  can  filtering  noise  near  the  edges,  thus  the  diffusion  tensor  which  controls 


the  diffusion  process  should  be  not  a  scalar  but  a  matrix  £> , 

anisotropic  diffusion  equation  is  introduced: 

dtu  =  div(D(Vua)-Vu) 

Here,  an  edge-enhancing  diffusion  model  is  adopted. 

Pi  =g( lViG  I) 


d\\ 

\d2i 


d  \ 

12  ,  and  nonlinear 

d22/, 


(3) 


(p2  -  0.2 


(4) 
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where  (px  controls  the  diffusion  in  the  gradient  direction,  cp2  controls  the  smoothing  in 
the  coherence  direction.  The  value  of  (p2  lies  in  the  interval  [0,  1], 


2.3  Diffusivity  Function  g 


Diffusivity  function  g  is  a  non-increasing  function,  its  value  lies  in  the  interval  [0,1]. The 
function  g (s)  plays  a  role  of  a  ‘fuzzy  detector’  of  the  presence  of  an  edge  at  a  particular 
position. 

Edge-enhancing  diffusion  can  adopt  the  diffusivity  function  proposed  by  Weickert: 

r  ~  \ 


g(s)  =  1  -  exp 


C 


{s!X)n 


(5) 


with  C=  -2.33667  and  m=4  here.  This  diffusivity  and  the  corresponding  flux  function  is  drawn 
in  Figure  1. 


Figure  1.  The  Weickert  diflfusivity(left)  and  the  corresponding  flux  function  s  ■  g(s)  (right) 


2.4  Discretization  Scheme 

In  order  that  the  continuous  anisotropic  diffusion  equation  is  applied  in  the  discrete  image 
data,  an  effective  semi-implicit  discretization  scheme,  e.g.  the  AOS  scheme,  can  be  employed. 

1  m 

uk+i  =  —  V  (/  -  m  tA,  (u  k ))  1  u  k  (6) 

m  ,= i 

with  1  is  the  direction  index,  1=1,..., m0 

The  AOS  scheme  will  separate  the  2D  diffusion  into  several  one-dimensional  diffusion 
process  along  chosen  directions,  it  is  absolutely  stable  for  any  discretization  time  step,  and 
contain  simple  and  effective  calculation  in  an  iterative  process. 

3  EXPERIMENT  RESULTS  AND  ANALYSIS 

3.1  Evaluation  Measures 
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For  objectively  evaluating  the  effectivity  of  the  nonlinear  diffusion  filter,  several  measure 
is  chosen  to  evaluate  the  despeckling  performance^  ljthe  capacity  of  despeckling  in  uniform 
areas;(2)the  capacity  of  preserving  edges  information;(3)the  capacity  of  keeping  the  average 
value  of  the  image. 

So  the  parameter  ENL  is  employed  to  measure  the  smoothing  perfonnance  of  the  speckle 
filter: 

ENL  =  n]  /  a\  (7) 

where  //  and  cr  represent  the  mean  and  the  variance  of  the  filtered  image  respectively.  While 

the  parameter  EPI  is  employed  to  measure  the  capacity  of  preserving  the  edge  of  the  filtering 
algorithm: 

EPI  =  Yj\  Ps  ~  Psn  \'Yj\Po-  Pon  I  (8) 

with  ps  is  the  pixel  value  of  the  filtered  image,  psn  is  the  neighborhood  pixel  value  of  ps,  p„  is  the 
pixel  value  of  the  original  image,  pon  is  the  neighborhood  pixel  value  of  p0. 

The  PM  represent  the  capacity  of  keeping  the  mean  value  of  the  image,  it  is  defined  as  the 
ratio  of  the  mean  value  of  the  original  image  and  the  filtered  image. 

3.2  Experiment  Analysis 

An  ERS-2  image  is  selected  for  the  despeckling  experiments.  The  experiment  test  and 
contrast  between  the  proposed  nonlinear  diffusion  approach  and  several  common  adaptive  local 
statistical  filter,  like  frost,  lee,  kuan  etc,  is  done.  Figure2-Figure7  show  the  filtering  results.  The 
parameter  used  in  this  diffusion  filter  is  shown  as  follow:  a  =  1.0,  r  =  1  ,T=2.  Table  1  give  the 
analysis  result  of  the  experiment. 

Firstly  in  terms  of  vision  effect,  the  performance  of  the  nonlinear  anisotropic  diffusion  filter 
in  SAR  image  despeckling  is  good,  it  not  only  effectively  reduce  noise,  but  also  preserve  the 
edges  information  of  the  image.  Seeing  from  its  filtered  image,  for  the  river  and  the 
homogenous  area  lying  the  river  side,  a  good  performance  is  shown,  and  the  texture  of  the 
mountains  appear  a  better  interpretation  perfonnance  by  eyes  than  the  original  image. 

From  table  1,  it  is  clear  that  for  any  evaluation  measure,  the  despeckling  algorithm  based 
on  nonlinear  diffusion  is  optimal,  others  algorithm  also  can  better  keep  the  mean  value  of  the 
image,  and  reduce  the  variance  of  the  image,  but  obtaining  a  good  perfonnance  both  in 
smoothing  noise  and  preserving  edge  appear  difficult.  However  nonlinear  diffusion  filtering 
algorithm  not  only  gain  a  good  balance  between  reducing  noise  and  preserving  edges,  but  also 
obtain  a  better  perfonnance  than  the  traditional  local  statistical  filter.  Furthennore,  because  of 
the  physical  essence  of  just  transporting  mass  of  the  diffusion,  its  capacity  of  keeping  the  mean 
value  of  the  image  is  so  strong,  the  mean  value  of  the  filtered  image  compared  with  one  of  the 
original  image  nearly  constant  according  to  table  1,  which  is  of  important  meaning  for  the 
further  image  interpretation.  Thus,  considering  the  evaluation  measure  all  together,  the 
nonlinear  diffusion  filter  can  effectively  reduce  the  speckle  in  the  SAR  images,  while  preserve 
the  important  detail  infonnation. 
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Table  1  The  contrast  of  several  despeckling  algorithm 


Mean 

Variance 

PM 

ENL 

EPI 

Original  image 

39.757342 

51.38114 

1.000000 

0.598725 

1.000000 

Kuan  Filter 

39.910042 

30.49218 

1.003841 

1.713119 

0.284019 

En-Lee  Filter 

39.567267 

31.527652 

0.995219 

1.575030 

0.313529 

En-Frost  Filter 

39.807742 

30.733126 

1.001268 

1.677728 

0.290327 

Gamma  MAP 

Filter 

38.095198 

30.769653 

0.958193 

1.532835 

0.296259 

AOS  ND  Filter 

39.757960 

30.298864 

1.000016 

1.722851 

0.330310 

Figure  2  Original  noisy  image 


Figure  3  Result  of  En-frost  filter 


Figure  4  Result  of  En-lee  filter 


Figure  5  Result  of  Kuan  filter 
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Figure  6  Result  of  Gamma  filter  Figure  7  Result  of  our  method 


4  CONCLUSION 

Speckle  is  the  intrinsic  characteristics  of  the  SAR  images  and  significantly  effects  the  remote 
sensing  application  of  the  SAR  images.  The  present  main  problem  of  the  traditional  adaptive 
filter  based  on  local  statistical  characteristics  lies  in  how  to  gain  a  better  balance  between 
preserving  texture  detail  and  despeckling.  So,  according  to  the  trait  of  nonlinear  diffusion  filter, 
the  attempt  to  introduce  in  SAR  image  despeckling  is  done,  and  by  experiment  test  and  contrast, 
the  proposed  approach  achieve  a  better  striking  experiment  perfonnance  than  the  traditional 
adaptive  local  statistical  characteristics  filter. 
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